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Air  Force  Technical  Applications  Center 
VELA  Seismological  Center 
Headquarters,  USAF 
Washington,  D.C.  20333 

Attention:  Lt.  John  Woods 

Subject:  2nd  Quarterly  Report  Covering  Period  1  December  1969 

to  31  March  1970 

Identification:  AFTAC  Project  VELA  T/0702/B/ASD 
Statistical  Discrimination 
ARPA  Order  624 
ARPA  Program  Code  9F10 

Contractor:  Texas  Instruments  Incorporated 

Contract:  F33657-70-C -031 1 

Effective  Date:  25  August  1969 

Amount:  $222,250 

Expiration  Date:  14  September  1970 

Project  Manager:  Stanley  J.  Laster 

Gentlemen: 

The  following  is  the  second  quarterly  report  of  statistical- 
discrimination  work  performed  under  Contract  F33657-70-C - 03 1 1.  The 
work  performed  and  plans  are  described  by  task,  with  more  details  being 
provided  in  the  attachment. 

I.  INTRODUCTION 

Research  under  the  current  contract  consists  of  the  following 

ta  sks: 

(1)  Continuing  use  of  the  Seismic-Event  Classification 
Software  Package  (SECSP)  at  the  Seismic-Array 
Analysis  Center  in  Washington 

(2)  A  literature  search  to  determine  the  observed  dif¬ 
ferences  between  the  seismic  signatures  of  nuclear 
and  chemical  explosions 


MAILING  ADDRESS  POST  OFFiCE  BOX  5621  • 


DALLAS.  TEXAS  75222  •  238-2011  •  CABLE  TEXINS 


(3)  Comparison  of  the  discrimination  capabilities  of 
the  Alaskan  Long -Period  Array  (ALPA)  and  the 
Large-Aperture  Seismic  Array  (LASA)  in  Montana, 
using  recordings  of  the  same  events  when  possible 
and  comparing  results  with  those  published  for  the 
Norwegian  Seismic  Array  (NORSAR) 

(4)  Processing  of  previously  recorded  events  from 
LASA  to  determine  the  optimum  procedure  for 
using  LASA  for  event  discrimination  and  study  of 
various  techniques  for  improving  long -period 
discriminants 

II.  WORK  SUMMARY  AND  PLANS  BY  TASKS 

(1)  TI  is  continuing  to  build  additional  ensembles  for  various  geographical 
areas  and  types  of  events  (e.  g. ,  the  Kurile  Islands,  deep  events, 
Novaya  Zemlya  explosions).  The  efficiency  of  various  combinations 
of  short-period  discriminants  is  being  tested,  and  use  of  combined 
long-  and  short-period  discriminants  is  just  beginning.  A  group  of 

29  new  short-period  events  obtained  from  SAAC  (three  of  them  also 
having  long-period  data)  should  provide  a  good  test  of  the  various 
"learning  patterns.  " 

The  contract  monitor  has  provided  a  list  of  30  events  to  be  run  through 
SEC SP.  All  of  these  appear  to  be  available  on  library  tape  at  SAAC, 
but  only  six  have  been  run  through  the  event  processor.  A  request  has 
been  made  to  have  the  other  events  processed,  so  the  running  of 
SEC  SP  on  these  events  will  be  delayed  until  this  has  been  done. 

(2)  The  literature  search  comparing  chemical  and  nuclear  explosions  is 
proceeding  slowly.  The  primary  interest  is  in  differences  observable 
at  teleseismic  distances;  there  have  been  very  few  such  observances 
of  chemical  events.  The  final  results  of  this  study  will  be  reported 
later  in  a  special  report,  so  this  topic  is  not  covered  in  the  attach¬ 
ment  to  this  report. 

(3)  Programs  have  been  written  to  measure  high- resolution  spectral 
peaks  in  short  time  gates  of  a  long-period  signal;  this  allows  group- 
velocity  measurements  for  each  peak.  The  higher  frequency  peaks 
are  presumably  associated  with  higher-order  surface  modes.  Hope- 
iuily,  such  energy  will  be  useful  for  discrimination  purposes.  Theo¬ 
retical  seismograms  generated  under  a  previous  contract  will  be 
used  to  test  further  the  method  of  isolating  higher-mode  energy. 
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(4)  The  matched- Altering  study  is  almost  complete,  and  chirp  filters 
of  appropriate  length  have  been  found  to  work  almost  as  well  as 
master  events  for  some  continental  areas.  In  addition,  simple 
measurements  of  group  velocity  give  a  reasonable  prediction  of 
how  effective  matched  filtering  is  for  a  given  event. 

Long-period  spectra  for  the  various  events  are  being  computed 
for  use  in  the  design  of  discriminants.  The  chief  need  now  is  for 
more  long-period  data. 

Action  by  AFTAC 

None 

i  inancial  Status 


Financial  status  will  be  presented  in  the  next  AMSR  to  be 
submitted  in  the  latter  part  of  April. 


Veiy  truly  yours, 


Stanley  J.  Laster,  Program  Manager 
Texas  Instruments  Incorporated 


SJ  L:ms 
Attachment 
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A.  SEC SP  PROCESSING 

The  Seismic-Event  Classification  Software  Package  (SECSP) 
is  a  set  of  programs  used  in  classifying  an  unknown  seismic  event  as  either 
an  earthquake  or  an  explosion.  Its  ultimate  purpose  is  to  detect  clandestine 
underground  nuclear  explosions. 

The  package  consists  of  two  phases  —  learning  and  classification. 
The  learning  phase  uses  statistics  from  an  ensemble  of  known  earthquakes  and 
explosions  (a  "learning  pattern").  Currently,  the  package  uses  statistics  from 
short-period  traces  only,  but  long-period  statistics  will  also  be  used  in  the 
near  future. 


The  learning  phase  develops  transformation  matrices  and  vec¬ 
tors  that  will  accomplish  certain  objectives.  Four  independent  pattern-recog¬ 
nition  methods  are  now  in  use: 

•  Class  clustering,  which  clusters  events  of 
like  class 

•  Class  separation,  which  separates  events  of 
unlike  class 

•  Likelihood  ratio  test,  which  estimates  the  prob¬ 
ability  that  the  given  event  belongs  in  class  1  or 
2.  (A  simplified  version,  linear  discriminant, 
assumes  the  same  covariance  matrix  for  each 
class  and  makes  a  decision  based  on  the  distance 
from  the  two  class  means.) 

•  Adaptive  algorithm,  designed  to  produce  a  vector 
whose  dot  product  with  the  event  statistics  is 
negative  for  one  class  and  positive  for  the  other 


A  - 1 


services  group 


The  statistics  now  being  used  are  P30  mean  square,  autocor¬ 
relation  mean  square,  envelope  difference,  and  spectral  ratio.  Bodywave 
magnitude,  dominant  period,  and  signal-to-noise  (S/N)  ratio  are  statistics 
to  be  tried  in  the  future.  Plan'*  also  call  for  testing  new  learning  patterns 
by  using  various  combinations  of  the  statistics.  The  various  statistics  have 
been  studied  by  means  of  2 -dimensional  plots.  The  P30  mean  square  and 
envelope  difference  (both  measures  of  complexity)  have  been  shown  to  be 
strongly  correlated  —  as  are  bodywave  magnitude  and  spectral  ratio.  Sur¬ 
prisingly  (because  of  the  similarity  in  method  of  computation),  the  P30  mean 
square  and  autocorrelation  mean  square  show  little  correlation.  In  addition, 
autocorrelation  mean  square  is  less  effective  than  P30  mean  square  as  a 
discriminant,  so  the  former  will  probably  be  dropped. 

The  initial  phases  of  the  processing  with  the  SECSP  programs 
used  a  learning  ensemble  of  16  earthquakes  and  10  explosions  assembled 
under  a  previous  identification  contract.  This  learning  ensemble  was  too 
small  to  give  an  accurate  statistical  representation  of  earthquakes  and  ex¬ 
plosions  and  did  not  represent  a  good  geographical  distribution  for  the  area 
of  interest.  All  except  one  of  the  explosions  were  from  a  single  area,  while  the 
majority  of  the  earthquakes  were  from  the  Kurile -Kamchatka  area.  In  such 
a  situation,  discrimination  could  arise  from  geographic  differences  rather 
than  from  actual  differences  between  explosions  and  earthquakes.  Thus,  a 
new  and  larger  ensemble  of  earthquakes  and  explosions  was  assembled  that 
hopefully  would  be  more  representative  of  the  area  under  study. 

The  new  ensemble  comprises  42  earthquakes  and  40  explosions 
selected  from  three  other  in-house  ensembles  in  addition  to  the  original  one. 
Due  to  the  limited  number  of  test  sites,  all  but  four  of  the  explosions  are  again 

from  the  same  area. 

The  earthquakes  were  selected  for  geographical  location;  they 
were  scattered  over  Russia,  China,  and  bodering  countries,  and  there  are  five 
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from  Northern  Africa.  Northern  Russia,  which  has  very  few  earthquakes,  is 
the  only  area  not  represented.  All  of  the  earthquakes  have  depths  of  <  70  km. 
This  learning  ensemble  will  be  used  by  subsequent  SECSP  runs,  although  16 
new  earthquakes  and  13  new  explosions  which  have  just  been  obtained  will  be 
added  as  soon  as  feasible. 

We  also  have  an  ensemble  of  35  earthquakes  deeper  than  90  km. 
Due  to  the  limited  number  of  in-house  deep  earthquakes,  all  deep  events  have 
been  taken,  regardless  of  location;  this  gives  a  worldwide  distribution  of 
deep  earthquakes.  Deep- earthquake  traces  resemble  explosion  traces  in 
many  respects  and,  therefore,  are  harder  to  distinguish  from  explosions. 
Working  with  the  deep  quakes  will  be  a  good  test  of  the  classification  program 
and  should  reveal  subtle  differences  among  the  three  event  classes. 

\  *  i  « 

V 

Power  spectra  are  computed  for  the  ensemble  of  shallow- focus 
events  (depth  <  70  km)  recorded  at  LASA.  The  ensemble  includes  39 
and  42  earthquakes.  The  average  spectrum  for  each  of  the  two  classes  of 
events  is  computed  and  shown  in  Figure  la.  The  spectra  are  normalized  to 
the  total  power  in  the  0.  85-  to  2.  80-Hz  band  before  averaging  inasmuch  as 
the  bodywave  magnitude  m^  ranges  from  3.  9  to  6.4  for  the  explosions  and 
from  4.  1  to  6.  1  for  the  earthquakes.  A  comparison  of  the  average 
spectra  is  indicative  of  the  effectiveness  of  the  spectral- ratio  discriminant  for 
distinguishing  between  earthquakes  and  nuclear  explosions;  e.g.  ,  the  average 
nuclear  explosion  generates  relatively  more  energy  in  the  1-3  cps  band  and 
relatively  less  in  the  .  5-1  cps  band  than  does  the  average  earthquake. 

Average  noise  spectra  computed  from  a  time  gate  preceding 
the  signal  for  the  two  classes  of  events  are  shown  in  Figure  lb.  The  spectra 
agree  closely  below  1.  4  Hz  —  but  above,  the  earthquake- ensemble  noise  is 
approximately  2  db  higher  than  the  explosion -ensemble  noise,  a  result  which 
would  tend  to  degrade  the  spectral-ratio  discriminant  for  low -magnitude 
explosions. 
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Next,  the  events  in  each  class  are  separated  by  magnitude 

(m,  <  5.  0  and  m,  >5.  0)  and  average  spectra  computed,  giving  the  results 
b  b 

shown  in  Figures  lc  and  Id.  For  this  ensemble  of  events,  the  effectiveness 

of  the  spectral-ratio  discriminant  is  largely  reduced  by  the  low  S/N  ratio  for 

events  with  m,  <5.  0. 

b 

B.  CONTINENTAL-PATH  STUDIES  -  DETECTION  OF  HIGHER-ORDER 
SURFACE-WAVE  MODES 

Preparation  of  the  necessary  computer  programs  has  received 
the  main  effort  in  the  task  to  detect  and  separate  any  Rayleigh  and  higher- 
order  shear-mode  energy  from  long-period  LASA  beams.  These  programs 
compute  and  plot  the  theoretical  dispersion  curves  from  layered  models,  the 
traveltirnes  for  the  various  wave  modes,  and  the  maximum  entropy  power 
spectra  of  the  time-partitioned  seismograms. 

The  theoretical  dispersion  curve  for  the  LASA  are  computed 

5jC 

fox  four  models:  Fort  Peck,  Til,  the  University  of  Wisconsin  Model  3  (UW3)f 
and  the  USGS  Model  3  (USGS3).  Dispersion  curves  for  these  models  are 
similar,  with  only  slight  differences  between  them  in  the  0.  0-  to  0.  1-Hz 
spectral  region.  Figures  2,  3,  4,  and  5  show  the  dispersion  curves  for  the 
Fort  Peck,  Til,  UW3,  and  USGS3  models.  From  the  dispersion  curves, 
theoretical  arrival  times  are  computed  for  the  Rayleigh  mode  and  the  first 
higher-mode  shear  waves  given  some  epicentral  distance.  While  it  is  not 
realistic  to  assume  that  the  LASA  model  applies  over  the  entire  great  circle 
route,  it  provides  a  starting  point  for  the  search  for  higher-order  modes. 

The  LASA  mod  ii  which  appears  to  give  the  best  results  in  preliminary 
analysis  will  be  used. 


Texas  Instruments  Incorporated,  1967:  Continuation  of  Basic  Research  in 
C rustal  Studies,  Final  Report,  ContractAF  49(638)- 1588,  AFOSR,  15  July. 
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Two  events  have  been  examined  for  higher -mode  content. 

Event  L-37,  a  teleseismic  explosion  (m^  =  6.  1)  has  a  signal*-power-to-noise- 
power  ratio  of  approximately  25  db.  Event  L-01,  an  earthquake  (m^  -  5) 

the  Greenland  Sea  with  an  epicentral  distance  of  approximately  5761  km,  has  an 
excellent  signal-power-to-noise-power  ratio  of  about  100  db.  The  vertical  com¬ 
ponent  of  the  long-period  seismograms  is  partitioned  into  blocks  containing  10G  or 

200  time  points  (100  or  200  sec  of  real  time).  The  first  block  was  taken  ap¬ 
proximately  300  sec  before  the  expected  first  arrival.  Each  successive  block 
of  points  overlaps  the  preceding  block  by  half  its  length  to  give  some  redun¬ 
dancy  to  the  computed  power  spectra.  The  maximum- entropy  power  spectrum 
is  then  computed  for  each  block  of  points  using  prediction- error  filter  lengths 
of  50  for  200-point  blocks  and  25  for  100-point  blocks.  The  frequency  and  power 
of  all  spectral  peaks  are  tabulated  and  assigned  an  arrival  time  equal  to  the 
middle  of  the  partition. 

Figure  6  shows  results  of  event  L-37  for  which  the  Til  model  is 
used.  The  spectral  peaks  picked  from  the  power  spectra  are  plotted  as  circles 
filled  by  an  amount  denoting  the  power  relative  to  the  largest  peak.  The  most 
powerful  peaks  tend  to  cluster  near  the  line  of  expected  arrival.  Such  good 
agreement  is  unexpected  because  of  the  assumptions  made  in  computing  the  ar¬ 
rival  time.  Inasmuch  as  L-37  is  a  subsurface  explosion,  little  or  no  second- 
order  shear-mode  energy  should  be  present.  What  appears  to  be  possible  sec¬ 
ond-order  energy  above  0.  05  Hz  and  between  1900-  and  2700- sec  traveltime 
may  not  be  significant  inasmuch  as  the  power  levels  are  consistent  with  the 
ambient  noise.  The  three  points  around  0.  04  Hz  after  3300- sec  traveltime  are 
definitely  above  ambient  noise.  Their  identification  must  await  further  analysis. 

The  results  for  event  L-01  are  shown  in  Figure  7.  The  Til  model 
is  used  for  the  traveltime  curve.  For  this  event,  the  mos*-  powerful  waves  tend  to 
arrive  earlier  than  expected  -  by  approximately  50  sec.  Above  0.  05  Hz,  the  dif¬ 
ference  increases  to  approximately  100  sec.  This  early  arrival  indicates  <-hat  the 
surface  layers  through  Eastern  Canada  have  a  higher  propagation  velocity  than 
do  the  LASA  models. 
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Figure  6.  Event  L-37  Power  Spectra  Showing  Spectral  Peaks 
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Figure  7.  Event  L-01  Power  Spectra  Showing  Spectral  Peaks 
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The  cluster  of  points  between  0.  06  and  0.  08  Hz  arriving  be¬ 
tween  1300  and  1700  sec  are  possibly  higher-mode  energy.  Having  energy 
levels  about  50  db  less  than  the  fundamental  mode  but  still  40  to  50  db  larger 
than  the  ambient  noise  level,  these  points  may  be  artifacts  generated  by  the 
analysis  technique  or  perhaps  arrivals  of  leaky  modes  not  normally  seen  but 
visible  here  because  of  the  large  S/N  ratio  of  this  event.  Both  possibilities 
are  being  investigated.  Theoretical  seismograms  of  known  content  and  dis¬ 
persion  generated  under  an  earlier  contract  will  be  tested  in  the  same  way  as 
the  events  to  see  if  these  peaks  are  real  or  anomalous. 

C.  LASA  PROCESSING 

1.  Group -Velocity  Analysis 

Long-period  surface  waves  from  a  group  of  explosions  and 
earthquakes  from  the  LASA  seismometer  arrays  in  Montana  were  analyzed, 
but  only  vertical  and  inline  components  were  used  because  the  transverse 
component  was  very  erratic.  Surface  waves  disperse  when  traveling  in  the 
earth  because  the  phase  velocity  of  propagation  of  surface  waves  depends  on 
wave  frequency;  the  original  impulse  spreads  or  disperses,  with  lower  fre¬ 
quencies  coming  in  before  higher  frequencies  in  the  frequency  band  of  interest. 

A  single  analysis  of  the  surface  waves  was  programmed.  The 
program,  operating  on  the  portion  of  the  trace  beginning  approximately  at  the 
predicted  Rayleigh-wave  onset,  uses  an  average  of  800  points  of  data  which 
are  optionally  bandpass -filtered  in  the  time  domain.  The  program  then  picks 
every  peak  and  every  trough  in  the  filtered  trace  and  assigns  it  an  order  num¬ 
ber  or  sequence  number.  At  each  point,  frequency  and  group  velocity  are  com¬ 
puted.  The  group  velocity  is  the  wavelet's  average  velocity  as  it  travels  from 
source  to  receiver.  Frequency  is  computed  by  numerically  differentiating  the 
order  numbe~  with  respect  to  arrival  time.  Finally,  grorp  velocity  vs  fre¬ 
quency  is  plowed  for  each  point.  The  study  used  26  tele  seismic  earthquakes 
and  11  tele  seismic  explosions  from  the  present  contract  ensemble.  All  explo¬ 
sions  are  from  a  single  area;  the  earthquakes  are  scattered  over  the  globe. 
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Analysis  results  vary.  The  explosions,  which  generally  have 
small  surface  waves,  show  the  expected  high  degree  of  scatter  and  only 
slight  trends  in  the  plotted  points.  This  result  is  due  largely  to  the  poor  S/N 
ratio  typical  of  explosions.  Figure  8  is  the  best  result  obtained  for  an  explo¬ 
sion.  Its  magnitude  is  5.  7.  In  addition,  earthquakes  with  low  S/N  ratios 
also  show  a  high  degree  of  scatter,  indicating  that  the  lower  the  S/N  ratio, 
the  poorer  the  results.  Figure  9  is  an  example  of  an  earthquake  with  low 
S/N  ratio. 

From  the  geographical  standpoint,  the  plots  reveal  that  events 
from  certain  regions  are  much  better  than  others.  Events  originating  in  he 

t 

Southwest  Pacific  show  very  good  results,  with  only  slight  scatter  (Figure  10  ). 
Better- than-average  results  are  also  found  for  events  in  Panama  and  South 
America  (Figure  11)  and  in  the  Greenland  Sea  (Figure  12).  At  the  opposite 
extreme  are  events  from  the  Kurile  Islands  area  (Figure  13);  events  from 
these  islands  and  from  neighboring  Hokkaido  and  Kamchatka  produce  plots 
with  a  h-gh  scattering  of  points.  The  crust  i.i  the  source  region  is  complicated, 
and  the  path  of  surface  waves  from  this  area  must  pass  through  the  complicated 
crust  structure  of  the  Aleutian  Islands  in  reaching  the  LASA  array;  this  is 
theorized  as  the  reason  for  the  poor  results. 

The  scatter  of  the  group- velocity  curves  should  indicate  roughly 
how  much  chirp- filtering  the  trace  will  enhance  the  peak  signal.  In  general, 
the  less  scatter  in  the  curve,  the  greater  the  signal  enhancement  so,  on  this 
basis,  we  predicted  that  chirp- filtering  a  Kurile  Islands  events  would  enhance 

>'fi 

the  signal  very  little;  workers  on  another  contract  found  this  to  be  the  case. 

The  group- velocity  curves  should  also  help  predict  the  average  crust  and  the 
average  dispersion  curve  over  the  path  of  the  surface  waves. 

# 

Texas  Instruments  Incorporated,  1970:  Seismic  Array  Processing  Techniques, 
Quarterly  Report  No.  2,  Contract  F33657-70-C - 0100,  AFTAC,  3  March. 

^Theoretical  curve  for  oceanic  path  taken  from  R.  Piermattei  and  A.  Nowroozi, 
1969,  Dispersion  of  Rayleigh  waves  for  purely  oceanic  paths  in  the  Pacific: 
Bull.,  Seismolcgical  Society  of  America,  v.  59,  p.  1913,  Fig.  9,  Oct. 
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Group  Velocity  Vs  Frequency  for  Earthquake  in  Easter 
China  as  Derived  from  Vertical  and  Inline  Traces 


Group  Velocity  Vs  Frequency  for  Earthquake  in  South 
Pacific  as  Derived  from  Vertical  and  Inline  Traces 


A-  18 


services  group 


Figure  11.  Group  Velocity  Vs  Frequency  for  Earthquake  in  Panama 
as  Derived  from  Vertical  and  Inline  Traces 


Group  Velocity  Vs  Frequency  for  Earthquake  in  Greenland 
Sea  as  Derived  from  Vertical  and  Inline  Traces 
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Figure  13.  Group  Velocity  Vs  Frequency  for  Earthquake  in  Kurile 
Islands  as  Derived  from.  Vertical  and  Inline  Traces 


2.  Matched  Filtering 

Application  of  the  M  vs  m  discriminant  to  shallow  teleseis- 

s  b 

mic  events  requires  a  good  estimate  of  the  long-period  Rayleigh  wave.  Gen¬ 
erally  speaking,  the  LASA  outputs  beamsteered  at  Rayleigh  velocity  show 
significant  waveform  degradation  due  to  poor  S/N  ratio  at  approximately 
4.  5  for  earthquakes  and  m^  <5.7  for  explosions. 

To  examine  certain  data-processing  techniques  for  enhancing 
S/N  ratio,  a  program  to  matched-filter  long-period  Rayleigh  waves  has  been 
written.  Matched  filtering  is  performed  by  crosscorrelating  a  matching  wave¬ 
form  (filter)  with  the  recorded  time  trace,  the  output  for  a  perfect  match  be¬ 
tween  the  filter  and  the  noise-free  time  trace  being  the  autocorrelation  func¬ 
tion  of  the  filter. 

Two  types  of  matched  filters  have  been  applied  to  an  ensemble 
of  explosions  (4.8  £  mfe  £  6.  1)  recorded  at  LASA.  The  first  is  a  master- 

event  (ME)  filter.  An  event  with  good  S/N  ratio  is  chosen  and  the  dispersed 
Rayleigh  waveform  partitioned  from  the  time  trace  for  use  as  the  matched 
filter.  The  program  truncates  the  filter  at  the  nearest  zero  crossing  and 
applies  an  8-point  cosine  taper  to  each  end;  in  addition,  a  zero-phase  band¬ 
pass  filter  may  be  designed  and  convolved  with  the  ME  filter  before  matched 
filtering.  This  filter  is  then  applied  to  other  events  with  low  S/N  ratio  occurr¬ 
ing  in  the  same  area.  The  second  type,  a  chirp  Liter,  is  generated  by  specify¬ 
ing  the  initial  frequency  fQ,  final  frequency  f  ,  and  time  duration  L,  where 


Both  techniques  have  been  previously  used  with  varying  success  by  many 
investigators. 
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Figure  14  plots  the  beamsteered  vertical  and  inline  traces  for 
the  events  used.  The  sample  rate  of  the  data  is  1  sec,  and  all  plots  presented 
in  this  section  are  scaled  by  setting  the  maximum  trace  deflection  to  1  in. 


The  vertical  Rayleigh- wave  component  of  event  L-27  (m^  =  6.  1) 
is  selected  for  the  ME  filter  (Figure  15)  because  it  exhibits  the  highest  S/N 
ratio  for  events  in  the  ensemble.  The  length  is  800  sec,  and  a  zero-phase 
high- cut  filter  at  0.  1  Hz  is  applied  to  the  ME  filter.  The  ME  filter  is  first 
applied  to  the  entire  time  trace  f*  om  which  it  is  taken  (L-37  vertical),  giving 
the  result  shown  in  Figure  15.  The  output  represents  the  best  S/N  ratio  im¬ 
provement  to  be  expected  from  the  ensemble  of  events.  The  S/N  ratio  is 
defined  here  as  the  maximum  zero-to-peak  amplitude  of  the  signal  divided 
by  the  rms  of  the  noise  preceding  the  signal.  The  ME  filter  is  also  applied 
to  the  inline  horizontal  trace  (Figure  15).  Ideally  this  output  would  exhibit  a 
90°  phase  shift  relative  to  the  vertical  trace  output.  An  S/N  ratio  improve¬ 
ment  of  approximately  2  over  the  input  trace  is  indicated.  The  beamsteered 
horizontal  traces  at  LASA  generally  show  a  much  lower  S/N  ratio  than  do  the 
corresponding  vertical  traces.  Results  of  applying  the  ME  filter  to  the  re¬ 
maining  events  (both  vertical  and  inline  traces)  in  the  ensemble  are  presented 
in  Figure  15;  the  events  are  ordered  by  decreasing  bodywave  magnitude. 

The  time  at  which  the  peak  in  the  matched- filter  output  is  expected  to  occur 
is  marked  for  each  event.  A  peak  in  the  outputs  of  the  matched- filter  vertical 
traces  is  obtained  for  events  down  to  m^  =  5.4  except  for  event  L-33. 

Next,  several  chirp  filters  have  been  designed  and  applied  to 
event  L-37  to  determine  the  best  match.  The  process  is  one  of  trial  and  error 
in  which  initial  and  final  frequencies  are  held  at  0.  025  Hz  and  0.  050  Hz,  re¬ 
spectively,  and  the  time  duration  of  the  filter  is  varied.  The  first  filter, 
which  is  800  sec  long,  results  in  the  output  shown  in  Figure  16;  in  this  case, 
the  chirp- filter  frequency  is  changing  too  slowly,  shifting  the  high  frequencies 
early  relative  to  the  predicted  peak.  Next,  a  600- sec  chirp  filter  is  applied, 
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the  low  frequencies  are  shifted  early,  indicating  that  the  chirp-filter  fre¬ 
quency  is  changing  too  fast.  Figure  16  shows  the  results  for  a  650-sec  filter  — 
the  best  match  obtained.  As  expected,  the  S/N  ratio  improvement  is  some¬ 
what  less  than,  that  obtained  from  the  ME  filter. 

In  Figure  17  are  the  15  individual  seismometer  channels  for 
event  L-37  and  the  matched- filter  results  using  the  650- sec  chirp  filter. 
Channels  F3,  D4,  Dl,  and  FI  show  poor  S/N  ratios  on  the  matched-filter  out¬ 
puts;  this  results  primarily  from  the  isolated  large  noise  excursions  seen  on 
the  input  traces.  The  noise  pulses  are  probably  nonseismic.  These  outputs 
indicate  four  channels  which,  if  used,  could  seriously  degrade  the  beamsteered 
output  and  subsequent  matched- filtering  results. 

For  the  remaining  11  channels,  the  peak  of  the  matched- filter 
output  is  timed  for  use  in  stacking  the  channels.  This  technique  yields  an 
optimum  beamsteered  output  as  opposed  to  one  using  theoretical  traveltimes, 
where  the  energy  is  assumed  to  propagate  along  the  great  circle  path  from 
epicenter  to  LASA.  The  11  channels  are  restacked  using  the  measured  delays 
and  filtered  using  the  650- sec  chirp.  The  output  is  shown  in  Figure  17.  The 
S/N  ratio  improves  slightly  over  that  for  the  chirp  result  shown  in  Figure  16, 
the  improvement  results  primarily  from  omitting  the  four  noisy  channels  from 
the  beamsteer.  The  measured  delays  are  close  to  the  theoretical  delays  for 
this  event,  the  maximum  difference  being  three  samples  (seconds).  The  11- 
channel  beamsteer  shows  a  4.  5-db  S/N  ratio  improvement,  and  the  matched- 
filter  beamsteer  shows  an  18.  5-db  S/N  ratio  improvement  over  the  raw  time 

traces. 

The  frequency  response  of  the  650-point  chirp  filter  is  shown 
in  Figure  18.  The  bandpass  characteristic  results  from  specification  of  the 
initial  and  final  frequencies  in  the  filter  calculation.  This  raised  the  question 
of  how  much  S/N  ratio  improvement  could  be  attributed  to  simple  bandpass 
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filtering  of  the  input  traces.  A  zero-phase  bandpass  filter  closely  approxi¬ 
mating  the  response  of  the  chirp  filter  was  designed  and  applied  to  the  L-37 
channels,  yielding  the  results  shown  in  Figure  18.  The  average  S/N  ratio 
improvement  is  2  db  over  the  raw  time  traces. 

One  additional  qu3stion  concerning  matched  filtering  is  how 

critical  is  the  design  of  the  matching  filter.  This  may  be  answered  by  seeing 

how  much  distance  can  be  allowed  between  the  epicenter  of  a  master  event 

and  the  epicenter  of  an  event  to  which  it  can  be  successfully  applied.  In  geo- 

* 

graphical  areas  where  matched -filtering  results  are  poor,  this  distance 
must  be  small  (<200  km),  but  this  need  not  be  true  in  other  areas. 

Shown  in  Figure  19  are  the  beamsteered  vertical  and  inline 
traces  for  earthquake  L-62  (m^  =  4.  3).  The  event  is  approximately  370  km 
farther  from  LASA  and  located  southeast  of  L-37.  The  ME  filter  from  L-37 
is  applied  to  event  L-62,  yielding  the  results  shown  in  Figure  19.  The  S/N 
ratio  improvement  is  about  6  db  over  the  raw  time  traces. 
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Texas  Instruments  Incorporated,  1970:  Seismic  Array  Processing  Techni¬ 
ques,  Quarterly  Report  No.  2,  Contract  F33657-70-C-0100,  AFTAC,  3  March. 
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